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Abstract 



A systematic analysis of large scale fluctuations in the low temperature 
pinned phase of a directed polymer in a random potential is described. These 
fluctuations come from rare regions with nearly degenerate "ground states". 
The probability distribution of their sizes is found to have a power law tail. 
The rare regions in the tail dominate much of the physics. The analysis 
presented here takes advantage of the mapping to the noisy-Burgers' equation. 
It complements a phenomenological description of glassy phases based on a 
scaling picture of droplet excitations and a recent variational approach with 
"broken replica symmetry". It is argued that the power law distribution of 
large thermally active excitations is a consequence of the continuous statistical 
"tilt" symmetry of the directed polymer, the breaking of which gives rise to the 
large active excitations in a manner analogous to the appearance of Goldstone 
modes in pure systems with a broken continuous symmetry. 
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I. INTRODUCTION 



The statistical mechanics of directed polymers in random media has attracted much 
attention in recent years Jl|-§[. This problem and related problems of higher dimensional 
manifolds are encountered in a variety of contexts, ranging from the fluctuations of domain 
walls in random magnets > to the dynamics of magnetic flux lines in dirty superconduc- 

tors . In addition, the randomly-pinned directed polymer is one of the simplest models 
that contain many of the essential features of strongly frustrated random systems such as 
spin-glasses [ID]. Understanding the behavior of the directed polymer is therefore important 
for developing intuition and testing theoretical ideas for more complicated random systems. 

Over the years, a variety of methods have been used to study random directed poly- 
mers. These include mapping @,^] to a hydrodynamic system: the noise-driven Burgers 
equation |lT|,12]|, the exact solution on a Cayley tree Migdal-Kadanoff approximate 



renormalization group calculations 14 , a Bethe ansatz solution in 1+1 dimensions using 



replicas [p^ , |lq1 , a gaussian variational ansatz in replica space ||17|| , and finally, renormal- 
ization group arguments and phenomenology |18| in the spirit of the droplet (or scaling) 



theory of spin-glasses [|19|]. There have also been substantial numerical simulations; some 
recent studies can be found in Refs. |I^ , pO| -[22[] . The qualitative phase diagram of the di- 
rected polymer is found to be quite simple: There is always a pinned phase dominated by 
disorder at low temperatures. For polymers in d + 1 dimensions, only the pinned phase 
exists for d < 2. But for d > 2, the polymer can undergo a continuous transition to a high- 
temperature phase where the disorder is irrelevant as has been proved rigorously ||23|| . Until 
recently, most of the efforts have been focused on characterizing the scaling properties of the 
polymer displacements and free energy fluctuations in the pinned phase. With the exception 
of 1+1 dimensions for which the scaling exponents can be computed exactly, systematic and 
analytic computations of the exponents at the zero temperature fixed point that controls 
the pinned phase have not been possible so far. 

Some of the other properties of the glassy pinned phase beyond the scaling exponents 
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were explored numerically by Zhang [p4| |, and more recently by Mezard [p0|| . These authors 



find very sensitive dependence of the polymer's low temperature configuration on the de- 
tails of the particular random medium. Such sensitivity is associated with rare but singular 
dependence on the details of the random potential of the system. This type of behavior, in- 
cluding sensitivity to small temperature changes, has been predicted by Fisher and Huse WE 



by phenomenological scaling and renormalization group arguments. As argued in Ref. |T8 



(and supported by numerical simulations in |18| , |20|^| ) , the physics of the low temperature 



phase is dominated by large scale, low energy excitations of rare regions, analogous to the 
droplet excitations proposed for Ising spin glasses ||19|| . Although the picture developed there 



is physically appealing and relatively complete, the phenomenological approach of Ref. [18| 



does not provide a systematic or quantitative way of calculating the properties of the pinned 
phase. On the other hand, various uncontrolled approximations can provide quantitative 



information. In particular, the Migdal-Kadanoff calculations of Derrida and Griffiths |14 
could be used to study the properties predicted in Ref. [[Tj|] although this has not been car- 
ried out in detail. Mezard and Parisi [T7| have recently proposed a very different approach: 



a variational method in replica space which can be used to study various aspects of the 
pinned phase. However, the method is limited by the gaussian Ansatz used whose physical 
significance is not clear, and the analysis based on replicas is haunted by the usual problem 
of the interchange of the n — * and the thermodynamic limits. This is particularly prob- 
lematic because within the gaussian Ansatz, one finds that the solution selected requires 
broken replica symmetry, the physical interpretation of which (if any) is unclear. (Note that 
the correct scaling exponents were obtained in 1+1 dimensions by using the Bethe Ansatz 
in replica space without breaking replica symmetry ||15|| .) 

In this paper, we explore the properties of the pinned phase of random directed poly- 
mers by using a more conventional approach based on a field-theoretic description without 
replicas, and the statistical symmetries of the problem. We will see that the existence of 
scaling forms for long wavelength, low frequency correlators of the noisy-Burgers' equation 
implies, without significant additional assumptions, the existence of rare large scale, low 



energy "droplet" excitations with a power-law distribution of their sizes. Although the large 
thermally active excitations are very rare, they dominate many thermodynamic properties 
and average correlation and response functions, as well as causing large variations in the 
properties of macroscopic systems. We therefore see that many of the properties of the 



pinned phase predicted by Fisher and Huse [113 and found in other uncontrolled approxi- 
mations can be recovered from the existence of a fixed point which, in the hydrodynamic 
language, is rather conventional. A fundamental lesson from this is that the broken con- 
tinuous statistical symmetry of the pinned phase of the directed polymer gives rise, quite 
generally, to power law distribution of large rare, low energy excitations; these are the analog 
for broken continuous statistical symmetries of the Goldstone modes associated with true 
broken continuous symmetries ! 

This paper is presented in a somewhat pedagogical fashion. It is intended to introduce 
the ideas of rare fluctuations as well as providing an alternative perspective for those familiar 
with the scaling approach to directed polymers and spin glasses. The paper is organized 
as follows: In Section II, we define the directed polymer problem and some of the glassy 
properties of the pinned phase. We motivate the considerations of almost degenerate ground 
states that give rise to large scale, low energy excitations, and relate their statistics to the 
distribution functions of the end point of a polymer. The distribution functions are computed 
in Section III by using a free-energy functional (described in Appendices A and B) and by 
exploiting the statistical symmetries (Appendix C). The results are interpreted in terms 
of the large, rare fluctuations in the pinned phase, with a short discussion contained in 
Section IV. The results of this study are then compared with the variational approach of 
Ref. ]17| and the phenomenological approach of Ref. [f[8| . We conclude that the various 
approaches point to the same physics governing the excitations of the directed polymer at 
low-temperatures — the physics of large, rare fluctuations. 



4 



II. PROPERTIES OF THE PINNED PHASE 



A. The Model 

We consider a directed polymer of length t (which will later be convenient to consider 
as "time") embedded in d + 1 dimensions. Let the position vector (£(z),z) with £ £ $l d 
and < z < t describe the path of the polymer in the d transverse dimensions. Then the 
statistical mechanics of this polymer is determined by the Hamiltonian 



Jo 



t 

dz 



(2.1) 



where k is the line tension, and r] is a quenched random potential of the medium through 
which the polymer passes. The random potential can be taken to be uncorrelated and 
gaussian distributed, with rj = and 



r]{x, z)rj(x', z') = 2D5 d {x - x')5{z - z 1 ) . (2.2) 

Here, the overbar denotes averages over rj, and D characterizes the strength of the random 
potential. Note that a cutoff on the spatial delta function is necessary for d > 2. 

In order to make our subsequent discussions precise, it is useful to fix one end of the 
polymer at an arbitrary point, say at = x. This is implemented by introducing the 
one-point restricted partition function, 

Z(x , t) = J V$ 5 d (i(t) - x) e-^VT } (2 3) 

where the arguments of Z give the coordinate of the fixed end point. Unless otherwise 
indicated, thermal averages will be performed using this restricted partition function, and 
will be denoted by (...). 

The partition function Z(x,t) is, of course, a random variable because of its dependence 
on the realization of the random potential rj. For a particular realization of r), the Hamil- 
tonian has no symmetries. However, because of the translational invariance in both x and 
t of the statistics of 77, the distribution of Ti. has statistical translational symmetry in x and 



t. Thus the statistics of Z and other properties will also be translationally invariant. As we 
shall see, this symmetry and the closely related "tilt" symmetry — which corresponds to a 
weakly t-dependent translation in x — have dramatic consequences. 

One of the simplest way to characterize the configuration of the polymer is the one-point 
distribution function: 

P 1 (y,to\x,t) = (5 d (i(t )-y)). (2.4) 

This function is the conditional probability of finding the segment, £(to), of the polymer at 
position (y, to) given that the fixed end is at (x,t). (We use "|" to separate the positions 
of the free | fixed positions of the polymer as in conditional probability.) In the absence of 
disorder, this distribution is easily calculated, yielding the usual random-walk result 



d/2 



G<°>(f -y,t-t ) = P?\y, t Q \x, t) = _ j exp 



re (x-y) 



2T t-t 



o 



(2.5) 



where the superscript (0) denotes the pure system with D = 0. The second moment of this 
distribution gives the mean transverse displacement, 

(|X(t-t )| 2 ) (0) ee (\i(t) -£(* )| 2 ) (0) = J d d rr-*G^(r,t-t ) = ™(t-to). (2.6) 

The existence of the random potential tends to increase the transverse wandering of the 
polymer, as it tries to take advantage of favorable regions of the random potential at low 
temperatures. After averaging over the disorder, translational symmetry in space is restored, 
and we have P\(y, to\x, t) = G t (x — y,t — t ). Note that a subscript t is used to indicate the 
explicit dependence of G on the polymer length; we have not put this subscript on P% since 
it already has the explicit t dependence. Numerically, it is found that 

(|A t (r)|2) ee J d d rr* G t (r, r) = r 2 <f(r/t), (2.7) 

with f(s) being a smooth scaling function which is finite throughout the range < s < 1, and 
the wandering (or roughness) exponent is 1/2 < ( < 3/4 depending on the dimensionality 



d ||25|| . In 1+1 dimensions, the exponent £ = 2/3 has been obtained exactly by a number 



of methods P,[15|. The power law dependence of the transverse displacement reflects the 
lack of intrinsic scales. This is possible because of the statistical translational symmetry, 
£ — > £ + const of the Hamiltonian ( |2 . 1| ) as discussed above. In dimensions d > 2, there 
is a high temperature or weak randomness phase in which ( = 1/2, in addition to the low 
temperature pinned phase on which we shall focus. It has been conjectured |H| that ( > 1/2 
for all finite d in the pinned phase although this is controversial. 

The detailed form of the disorder-averaged end point function Gt is much harder to obtain 
than the scaling exponent, even numerically pl|,22| . It is expected to have the general form 



G t (f,r)^r- d % T (r/r<) (2.8) 

by simple scaling and normalization requirements, with the scaling function g depending 
only weakly ont/r but decaying rapidly for r = \r\ ^> r^. Two limits of particular interests 
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G(f,r)= lim G t (f,T)«T-Scx>(r/T C ), (2.9) 

t/T— tOO 

which describes the one-point distribution of a semi-infinite polymer, and 

G{r,t) = G t {r,t) « r d %(r/t<), (2.10) 

which describes the distribution of position of the free end of a finite polymer. As ex- 
plained in the Appendices, analytical studies of the problem is much simplified in the limit 
t/r — > oo. For instance, G{r,r) is given to a good approximation by a self-consistent inte- 
gral equation in 1+1 dimensions (See Appendix B and Ref. ||26|| ). Numerical solution of the 
integral equation |26| yields the form Eq. (|2.9| ), with goo{s) well approximated by a gaussian, 



although the precise shape of the tail is not known. However, in numerical simulations, it 



is most convenient to study the end-point distribution G(f,t). It is found pl| , p2| that the 
distribution function obeys the scaling form Eq. ( |2.10| ). A comparison between goc{s) and 
gi(s) shows that the scaling function at the two limits are quite similar: both are sharply 
decreasing functions whose widths are of order unity. Thus we see that the semi-infinite 



polymer result G(f, r) gives the correct scaling behavior and the qualitative form of the 
scaling function for the end-point distribution of a finite polymer. We shall use this approx- 
imation in the following sections, where we will compute various distribution functions for 
a semi-infinite polymer (i.e., for r <C t) and then apply them to the end point of a finite 
polymer with r = t. 

B. Ground States 

The description of the directed polymer given so far is rather conventional. The form of 
G t in Eq. Q2.8|) describes a generic polymer. To appreciate the glassy nature of the pinned 
phase, it is necessary to go beyond the description in terms of the mean Pi. We first note 
that an exponent £ > 1/2 in the pinned phase implies that the energy scales involved in the 
pinned phase are large. From the first term in ( |2.1| ), we see that a displacement of order 
of the free end costs a minimum elastic energy of order t 2C -~ l which grows for long polymers. 
Growth of the characteristic energy scale for order parameter variation with length scale 
in conventional systems implies that the system is in an ordered phase governed by a zero 
temperature fixed point. Although there is no "order parameter" in the directed polymer, the 
displacement £ (z) play a similar role and growth of the energy scale for variations of £ with 
length scale — here t — implies that the pinned phase is controlled by a zero temperature 
renormalization group fixed point whose properties control the scaling of various quantities 
such as G t in Eq. (|2.8j ). By analogy with the ordered phase in conventional statistical 
mechanical systems, we expect that the large scale properties of the pinned phase can be 
described in terms of a "ground (or equilibrium) state" or "states" , and fluctuations about 
or between these "states". Thus the configuration of the polymer selected by "thermal" 
averaging should be the equilibrium "state" that optimizes the total free energy for a given 
realization of the random potential. At zero temperature for a finite polymer with one point 
fixed at (x, t), there will be a unique preferred path (see Fig. 1(a)), £*(z)> which is the ground 
state, i.e., the state with the lowest energy. Here, the term "state" refers to an optimal path 



S 



starting from the fixed end at (x, t). A well-defined thermodynamic limit exists for the state 
of a semi-infinite polymer if the thermal mean position (£(2)/ an d all other properties of 
the polymer at fixed finite z/t tend to a unique limit for a specific sample as t —>■ oo. The 



conjecture that this holds for almost all samples was made and supported in Ref. [{IS ]. The 
above definition of state only makes sense at T = after providing short distance cutoffs 
and a T in the £ and z directions respectively. At small but finite temperatures, thermal 
fluctuations wash out the effect of disorder at short length scales, and the polymer does not 
"feel" the random potential until scales a%(T) and a T (T) ~ a^ 2 {T) [fl8l| . For example, we 
have in d — 1 

a € (T) = T 3 /kD. (2.11) 

We can then study the "states" of the polymer coarse-grained on the scales a^(T) and a T (T). 
A natural conjecture is that at finite temperature the equilibrium state is still unique. What 
does this mean ? If the one-point distribution function, Pi(y, to\x, t) for a typical sample 
has the general behavior sketched in Fig. 1(b) (solid line), i.e., a sharply-peaked function 
of y centered about some y* = £*(£o)> with a width of order a^(T) for r = t — t ^> a T (T), 
then this, together with similar behavior for all < z < t, implies that there is a well- 
defined "state" £*(z) at long length scales at finite temperature. As we shall see, this simple 
picture is roughly correct but with subtle modifications of the meaning of "typical" which 
crucially affect the physics. Note that the equilibrium state is the minimum of the coarse- 
grained Hamiltonian which includes effects of the entropy of small scale fluctuations. Thus, 



in general, as shown in Ref. ||18|| , the equilibrium states for different temperatures of the 

same sample will be very different on long length scales. 

The peak in Pi (solid line in Fig. 1(b)) should be located within a transverse distance 
of x for most samples, since the disorder- average Pi, sketched as the dashed curve in 

Fig. 1(b), has the characteristic scale \y — x\ ~ r 1 *. The sharpness of P\ for a typical sample 

can be revealed by probing disorder moments of P\ for different end points. For instance, 

from the joint two-point correlation function, 



Qt(yi - x,y 2 - x,t - t ) = Pi(y 1 ,t \x,t)P 1 (y 2 ,t \x,t), (2.12) 
we can obtain the mean square of the (thermally averaged) displacement, 



(X t {r)) = J d d y x d d y 2 {y x - x) ■ {y 2 - x) Q t {& -x,y 2 - x, r), (2.13) 
which is expected to behave as 



Mr)) 2 ^r^f(r/t), (2.14) 

just like (^\X(t)\ 2 ^ in Eq. (|2.7| ) in the pinned phase of the directed polymer but with a 
different scaling function /'. From Eqs. ( |2.7|) and Q2.14j ), we obtain the mean square of the 
thermal fluctuations about (£(t)), 



C T (T,t)= (ji(t - r) - (i(t - r))f) = (|* t (r)|») - (x t (r)) 2 

= \^ d d y x d d y 2 \y x -y 2 ^Q t {y x -x,y 2 -x,r). (2.15) 

One would hope that this quantity Ct(t, t) would characterize the typical width of the one 

point distribution Pi (solid line of Fig. 1(b)). From the above discussion, the conjectured 

r— 1 1/2 

uniqueness of the equilibrium states at finite T would suggest that Ct(t, t) would be 
the diameter of the "tube" in which the polymer fluctuates, with Cxir^t) not growing with 
t or t. We shall see however, that in fact Ct{t, t) ~ r in both the pinned and the high 
temperature phases. How is this consistent with the existence of an equilibrium state of the 



problem in the pinned phase ? It has been argued |18| that for typical samples (in fact in 
almost all samples for large r) Ct{t, t) is of order unity, but rare samples have sufficiently 
large Ct{t, t) that they dominate the sample (or disorder) average. In this paper, we will 
see that this conclusion arises in a very natural way and can be characterized analytically 
by a power law tail in the distribution of various quantities that are related to Ct{t, t). 

In conventional critical phenomena, it is often sufficient to characterize a system by 
characterizing a few moments of the fluctuations. However, for many random systems, such 
averages do not give an adequate account of what actually goes on in a typical sample, 
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since some properties are not "self-averaging". The pinned phase of the directed polymer 



exhibits exactly this type of behavior p!8| , |20|j2^] and thus cannot be well characterized by the 
knowledge of a few moments. For instance, we shall see that a small fraction of samples have 
nearly degenerate ground "states", an d £,b(. z )> with segments (say the free ends £„(0) = 
y* a and ££(0) = yl) located far away (i.e., \y* a — y£\ ^> a^(T)), as shown in Fig. 2(a). Such 
states would be manifested in Pi as widely separated peaks (Fig. 2(b)) at low temperatures. 
The occurrence of such degenerate "states" may be rare. However, if they survive in the 
thermodynamic limit, with \y* a — yl\ — > oo as t — > oo, then they can still have important 
physical consequences, because these are just the type of states that give rise to large scale 
fluctuations at low temperatures. Such almost degenerate states are the low dimensional 
analog of the "droplets" [0,0 conjectured to govern much of the low temperature properties 
of spin glasses. The investigation of these large fluctuations is the focal point of this study. 
We are particularly interested in obtaining the statistics of the almost degenerate states and 
understanding the effects of these states on physical observables. We shall give a detailed 
account of these effects in the following sections. 



C. Distribution Functions 

There are a number of ways one can characterize the rare fluctuations. To find the 
relative abundance of samples which behave as Figs. 1 and 2, we can, for example, consider 
the probability that Pi has the structure shown in Fig. 2, with two peaks separated by a 
displacement A = |A|at,2 = £ — r, i.e., at a distance r from the fixed end. It is useful 
to consider the larger peak to correspond to the optimal path and the other to correspond 
to an excitation with displacement A away from the optimal path. For the second peak to 
have a reasonable amplitude, the excitation free energy must be of order T or less, i.e., it is 
a thermally active excitation, which is analogous to an active "droplet" introduced for spin 
glasses. The probability of an active droplet excitation of displacement A at a distance r 
from the fixed end is obtained from the distribution Q t as |37| 
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W t (A,r) = jd d y 1 d d y 2 S d (y 1 -y 2 -A)Q t (y 1 -x,y 2 -x,r). (2.16) 

Like the one-point distribution function G t discussed in Sec. II. A, it will be convenient to 
consider 

W(A,r) = lim W t (A,r) (2.17) 

t/r— >oo 

which, for large r, is the probability distribution of active droplets in the bulk of a semi- 
infinite polymer. Note that the somewhat different distribution of active droplets at the free 
end (z = 0) 

W(A,t) = W t (A,t) (2.18) 

is however the most readily measurable quantity numerically. 

In the absence of disorder, the "bare" distribution is just a gaussian, 

W (0, (A,r) = ( i ^)" / %-* il (2.19) 

since 

q(o) = p (o) . p (o) independent of t [ Note that Wt (£, r) is the probability that the 

— * 

distribution Pi has weight at two points separated by A. Only if the samples in which 
this occurs behave as in Fig. 2 will the interpretation of this as the probability of a well- 
defined active excitation really be useful. In the absence of randomness, this interpretation 
is clearly incorrect.] With disorder, Q t becomes highly nontrivial, and it is not easy in 
general to compute even the first few moments of W t . The large- A tail of the distribution, 
which governs the large scale fluctuations at low temperature, is difficult to obtain both 
analytically and numerically. In Section III, we introduce a free-energy functional which 

— * 

enables us to compute the asymptotic form of W(A, r) explicitly. We will show that the 
tail of W(A, t) [and hence also the tail of W(A, t) if we ignore the dependence on t/r] has 
a power law form. 

One of the physical consequences of the nearly degenerate ground states is the behavior of 
the mean square thermal fluctuations of a given segment of the polymer, CV(t, t) = (A 2 ) /2, 
which has an average value 
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C t (t, t) = \& = \j d d AA 2 W t (A, r). (2.20) 

This quantity is related to the susceptibility of the polymer to a tilt, i.e., the response to a 
term added to the Hamiltonian ( [2.1| ), 

H R [l V }=H[£, V }- fdzh-f z . (2.21) 



to 



Thermal fluctuations of £(£o) a t a distance t = t — t from the fixed end are simply propor- 
tional to the linear response [^] of the polymer by the fluctuation-susceptibility relation, 



x[v] = -S^Ut) - 6(tb))ji_o = ^C T (t - t ,t), (2.22) 

— * — # 

where ^ is a component of £, and the subscript ft, denotes thermal average taken with respect 
to the new partition function, 

Z(x, t-h)= j 6 d (C(t) - x) e- n rJ?. (2.23) 

Since the random potential r\ in the Hamiltonian Q2.1p does not single out any preferred direc- 
tion in the (£, z) space, a statistical tilt symmetry — related to the translational symmetry - 
exists which is recovered upon disorder average. As a result, one can prove straightforwardly 
that on average, the applied tilt h merely shifts the mean position of the polymer at to to 



£(to)y- = x — hr j k (see Appendix C and Ref. p9[). The disorder-averaged susceptibility 



is thus just 



which is exactly the same as the susceptibility of the pure system. 

If the statistical properties of the ground states are complicated as illustrated in Figs. 1(a) 
and 2(a), then the mean susceptibility x does not provide an adequate characterization. On 
the one hand, a sample such as the one shown in Fig. 1(a) will contribute very little to 
the average susceptibility at low temperatures since it is locked in a unique state, separated 
from the lowest excited state by a free energy difference ^> T. On the other hand, a sample 
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corresponding to Fig. 2(a) has nearly degenerate states. Thermal fluctuations will thus cause 
large scale "hopping" of the polymer from one state to the other. Thus the latter samples 
can give large contributions to the average susceptibility even though they occur rarely. The 
relative abundance of such samples is given by the susceptibility distribution function, 

D t (x,r) = 6(x-x[v})- (2-25) 

A first principles calculation of D t would require the knowledge of the full distribution of 
the function P 1; or at least all of its correlations, 

Pi(jfi,to\S,t) ■ ■ ■ Pi(y n ,t \x, t); 

this will not be attempted here. To obtain the qualitative behavior, we will instead assume 
that a sample with an active droplet of size A has a susceptibility of order \ ~ A 2 /(2Td). 
This assumption obviously will not be very good if there are many almost degenerate states 
that contribute equally to the susceptibility, but it should be reasonable if large scale de- 
generate states occur only rarely, as we will demonstrate is the actual case. Using this 
approximation, the susceptibility distribution can be simply obtained from the droplet dis- 
tribution as 

A(X, r) = j J^ Td) W t (A, r) u[ X /(A 2 /2Td)], (2.26) 

where u(s) is a sharply peaked scaling function whose precise shape depends on more knowl- 
edge of the distribution of Pi. 

D. Free Energy Variations 

The functions W t (A,r) and D t (x,T) introduced above describe the distribution of the 
polymer's equilibrium fluctuations and susceptibilities. But in reality (experimentally or 
numerically), equilibration of a glassy system is often very difficult ||30|| . The reason is that 
the nearly degenerate states that give rise to large scale, low energy droplet excitations are 
typically separated by large energy barriers, and therefore have extremely long relaxation 
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times. Knowledge of the distribution of the barriers is therefore crucial in understanding 
the dynamic properties of the glassy polymers. 

The details of the dynamics at low temperatures is very complex and beyond the scope of 
this paper. However, one can get a bound of the free energy barriers within the equilibrium 
theory. To do so, let us consider the nearly degenerate states shown in Fig. 2(a). To probe 
the free energy "landscape" between and y£, we would like to know F(x,t;y,0), the free 
energy of a polymer £(z) with both ends fixed at £(t) = x and £(0) = y, and then vary y in 
between y^ and yl (see Fig. 3). The maximum, F*, of F(x, t; y, 0) in this range of y gives a 
lower bound of the free energy barrier to move the end £(0) from y^ to y£, since the polymer 
must pass through all intermediate states with y* a < y < y£, i.e., in a range of width A. Note 
that this does not take into account the additional barriers the polymer may encounter at 
intermediate states (see Fig. 4). At this point, it is not clear whether or not the cumulative 
effects of such barriers will be much larger than 5F* = F* — F(x, t; y^, 0). 

Unfortunately, with two ends fixed, F is difficult to compute analytically. What can be 
obtained readily (see Appendices A and B) is the statistics of the free energy of a polymer 
of length t with only one end fixed, which is just 

F(x,t) = -T log Z(x,t). (2.27) 

As shown by the numerical studies of Kardar and Zhang [f|] , two polymers with end points 
fixed at a distance A < t** apart will merge a distance A 1 ^ from the fixed ends and coincide 
for the rest of the way (see Fig. 5). Since the free energy difference between two polymers 
with different end point positions will only arise from the section over which they differ, we 
expect the typical free energy difference, F(x, t; y, 0) — F(x, t; ife, 0), of an intermediate state 
in Fig. 3 with \y — y^\ ~ O(A), to scale the same way as the difference F(x + A, t) — F(x, t) 
for the two polymers in Fig. 5. The latter is characterized by the free energy correlation 
function 

C F (A, t) = [F(x + A,t)-F(x,t)\ 2 (2.28) 
which is conjectured to scale as 
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C F (A,t)~A 2a (2.29) 

for A <C ft, with the exponent known to be exactly a = 1/2 in 1+1 dimensions. Thus 
\J~C~p{K~t) provides the lower bound for the barrier of formation of an intermediate active 
droplet, which is the typical scale for F*. 

On the other hand, two polymers with end points fixed at A will not overlap at all 
and thus behave as if they are in independent samples. In this case, 

C F {A,t)~t 2d (2.30) 

with 9 being the exponent characterizing the sample-to-sample free energy variations. A 
simple scaling form, which we shall see arises naturally, links the two limits: 

C F (A,t) = A 2a c(A/t<) (2.31) 

yielding the scaling relation 

a = 6/C (2.32) 

We will later derive the scaling law 

9 = 2C - 1 (2.33) 

which arises simply from the naive contribution to the free energy from the elastic part of 
the Hamiltonian with displacement of order t^. 
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III. STATISTICS OF RARE FLUCTUATIONS 



In Section II, we defined a number of distribution functions and disorder averaged cor- 
relations functions which are useful in probing the glassy nature of the randomly-pinned 
directed polymer. To obtain these functions in a systematic way, and to uncover the inter- 
connections among them, we shall use field theoretic methods. The disorder average will be 
replaced by an average over a weighting functional, which describes the probability distri- 
bution of the free energy. This method is inspired by the Martin-Siggia-Rose dynamic field 
theory |31|] developed in the context of stochastic dynamics onto which the directed polymer 
can be mapped The explicit mapping to the noisy-Burgers equation is derived in Ap- 
pendix B. Note however, that the method may not be limited to the directed polymer and 
is described in Appendix A for an arbitrary dimensional manifold for which the mapping to 
stochastic dynamics cannot be performed. In this section, we shall use this field theoretic 
method to obtain the tails of the distribution functions introduced in Section II. C. We shall 
limit our discussions to the semi-infinite polymer (t — > oo) problem for which the analysis 
is the simplest. We then interpret the results to characterize the statistics of the large rare 
fluctuations. 



A. Derivation of the Distribution Functions 

1. Formalism 

To obtain the tails of the distribution functions W t (A,r) and D t (x, t) defined in Sec- 
tion II. C, we need (see Eq. ( p,16|) ) the disorder-averaged function Qtijji—x, y 2 — x, r) defined 
in Eq. (|2 . 1 2| ) . As shown in Appendix A, various distribution functions can be generated by 



adding a source term to the Hamiltonian, e.g., 

n^n+ /* dzl(i(z),z), (3.1) 

Jo 

and then differentiating the corresponding average free energy, F(x, t; J) with respect to J. 
The resulting expressions are simple for the semi-infinite polymer. For instance, 
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' )! ' U - ' J) - P x %to\x,t) = G(x-y,t- t ), (3.2) 



SJ(y,t ) 
and 

S 2 F(x , t; J) 



G2 i(x — y±, x — i/2) t — t ) = 

8J{y u tQ)5J{y2,t Q ) 



r (f(*o) " ft) 5 d (f(*o) - ft) ) - Q(Vi -x,y 2 -x,t- t ) 

In terms of the noisy-Burgers equation, the limit t —>■ 00 corresponds to the statistical 
steady state, and G and G 2 ,i are the linear and nonlinear response functions respectively (see 
Appendix B). It is instructive to consider the meaning of G 2j % in the limit of zero temperature 
(with an appropriate short distance cutoff on the random potential): G 2 ,\{x — ft, x — y 2 , r) 
is the limit of small e of 1/e times the probability that the optimal paths passing through 
yi and 1/2 a distance r from the fixed end at (x, t) both have energy within e of the ground 
state. Since triple degeneracies are unlikely, this is non-zero when the ground state is doubly 
degenerate and G 2 ,i is the "density of states" of degeneracies. In terms of G 2 ,i, the droplet 
distribution is 

W(A,t) = 5 d (A) + W'(A,r), (3.4) 

with 

W'(A, r)=T J d d y x d d y2b\y x - ft - A)G 2 ,i{x - ft, £ - ft, r). (3.5) 

The free-energy functional described in Appendix A allows us to express higher-order 
distribution functions such as G 2j i in terms of the one-point function G (see Eq. (|A16 )). For 
the directed polymer, the relevant result is 

G 2 ,i(x — y x ,x — y 2 ,t — t ) = — J d d x! ' d d y' 1 d d y 2 dt 'dt 1 dt 2 G(x — f ,r - t') 

ri, 2 (^ - ifi, t' - h; f - tf 2 , t' - t 2 )G(y 1 - jft, tt - t )G(y 2 - t 2 - t Q ), (3.6) 

where is a "vertex function" which is natural in the context of the noisy-Burgers equa- 
tion; it will be specified shortly. The expression can be represented diagrammatically as 



in Fig. 6. We see that the joint distribution G?2,i has been conveniently broken up as a 
convolution of a number of one-body distribution functions (the G's), with the branching 
process controlled by the vertex Ti^. 

For a (statistically) translationally invariant system, it is convenient (after disorder av- 
eraging) to work in Fourier space, with 



G(q, t)= J d d fG(f, r)e iq - r . (3.7) 
The Fourier transform, W, of W in Eq. (|3.5|) becomes 



W'{q, t-t ) = -T f dt' f dt x f dt 2 G{0, t - t' 

JO Jtn Jtn 



'0 J to J to 

f iM t'-h; -q, t! - t 2 )G{q, h - t )G(-q, t 2 - t ), (3.8) 



where r 12 is the Fourier transform of the vertex function r 12 . From the scaling form Q2.9| ) 
for G, we have 

G(q,r)^g(q^), (3.9) 

where q = \q\. Also, G(Q,t) = 1 by normalization of the probability distribution P\. There- 
fore we have the simple expression 

W'(q, r) = — T f dr' f d n f dr 2 f h2 (q, r' - n; -q, r' - T 2 )g(qr^g(q4). (3.10) 

J-oo JO JO 

To procede further, we need to know something about the vertex r 12 . We are especially 
interested in the small q behavior since that is what will control the tail of the distribution 
W(A, r). Note that the normalization of W requires that W'(q — » 0,r) — > 0, hence we 
must have T 12 (q —> 0) — > since g(0) = 1. Also, we recall that the mean susceptibility 
is x = T / K due to the statistical tilt symmetry (Appendix C), and x = Wd^- 2 from the 



fluctuation-susceptibility relation ( p.24|) . Thus 



_ 1 9T 

lim lim -VlW'(q, r) = -A 2 = — r. (3.11) 

T->oo q ->0 d K 



Since Vg-g = by symmetry, we must have from Eq. ( p. 11 
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f i, 2 (g, t' - n; -g, t' - t 2 ) = -g 2 7(g, t' - n; -g, r - r 2 ), (3.12) 

with 

dr' / dri / dr 2 7(5, r' - n; -g, r' - r 2 ) = r. (3.13) 
-00 Jo Jo 

The simplest form of r 12 satifying Eq. ( |3.13| ) is 

rPj(?y - n; -g>' - r 2 ) = ig 2 5(r' - n)5(r' - r 2 ). (3.14) 

This is actually the form of the "bare" vertex, i.e., it is the exact vertex in the absence of 
the random potential, as the bare distribution in Eq. (|2.19|) is readily recovered by 
substituting Eq. fl3~T4D in Eq. (|3~T0|) . In Appendix C, we use the statistical tilt symmetry to 
derive some Ward identities which ensure that the scaling behavior of the droplet distribution 
obtained from the full vertex r 12 will have the same form as that obtained by using the 
bare vertex fj° 2 . Using Eq. (|3.14j ) then, Eq. (|3.10| ) becomes 

2 T 



T 



K JO 



W'(q,T) w --q 2 I dr' ^(g(r') c ) " « --q^^w^T), (3.15) 



K 



where w(s) is a scaling function whose precise form depends on the actual forms of g and I\ 2 
(see Appendix C), but whose limits are simple, i.e., w(s — > 00) = const and w(s — > 0) — > s. 



2. Results 

The Fourier transform of W itself, W = 1 + W, has a power-law singularity (an inverted 
cusp) for small q in the limit r — > 00 as long as ( > 1/2. The singularity is only cutoff by 
r, with W ~ 1 — (T/n)q 2 T for g < r - ''. Inverse Fourier transforming VF, we obtain the final 
result 

W(A,t)k Ad+ l_ c _M ^/T C ) for A»a 5 (T), (3.16) 

with the scaling function wis) having the form w(0) = const and w rapidly decreasing for 
s ^> 1. This result suggests the following scaling form for the full distribution for a polymer 
of length t, 
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W t (A,r) 



(3.17) 



with w = Woo. Assuming that W t has a weak dependence on t/r as in G t , the above result 
leads to a power law distribution of active droplets at the free ends, 

1 



W(A,t) 



(3.18) 



which is cut off only by the finite length of the polymer. A similar result is obtained for the 
tail of the susceptibility distribution. From Eqs. ( |2.26| ) and ( |3.17D , we have, 

1 



D t (x,r) 



X 



2-(2C) 



-d t /r(x/r 
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for x ^ CLc/Td, 



(3.19) 



with d t / T being another scaling function which is qualitatively similar to w t / T . 

Note that Eqs. (|3.16|) through ( 3.19 ) should only hold in the scaling limit A 3> a^(T), 
where the one-point function has the scaling form (|3.9|) . To understand the form of crossover 
of W(A,t) from A < to A > a^, it is useful to consider the following approximate form 
of the one-point function, 



G(q, t) = e 2, 



: v{q)q 2 T 



(3.20) 



where V(q) = 1 + (qa^)^ 1_2 . This simple form of G extrapolates smoothly between the 
bare function G (q,r) = e~^ q2r for qa^ ^> 1 and the appropriate scaling form ( |3.9p for 
qa^ < 1. It is therefore a useful guide to the qualitative features of the crossover. (A better 
determination of the response function is given in Ref . f26| . ) Using the bare vertex Eq. (gl4|) 
and the approximate form Eq. (|3.20|) for G, the droplet distribution becomes 



W(q,t) = l- 


%) ^ 




(3.21) 


with the limiting forms 












1- v 


for 


q < T~t, 




W(q,r) = < 




for 


a^ 1 > q > r _c >, 


(3.22) 






for 


q > dg 1 . 
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Fourier transforming leads to a droplet distribution W(A,r) sketched in Fig. 7 (solid line), 
with the asymptotic scaling behavior for A > given by Eq. (|3.16|) , and a smooth behavior 



for A < a,£. It should be emphasized that the precise form of G and r 12 will change only the 
details of the crossover function connecting the scaling region to the regions with A < and 
A > t\ but not the qualitative features of the distribution. Assuming weak dependence of 
W t (A, r) on t/r, we expect the distribution W(A, t) for the end point to have the same form 
(solid line of Fig. 7) with yet another crossover function. Finally, the same considerations 
lead to a similar form for the distribution of the susceptibility, D t (x, T )- 

B. Glassy Properties of the Pinned Phase 

1. Droplet Excitations 

We now interpret the results obtained in the previous sections in terms of the structure of 
ground states and excitations of the randomly-pinned directed polymer. We first compare 
the result for W(A, r) in the limit of zero disorder. In this case, the microscopic cutoff 
length diverges, i.e., — > oo, and the distribution W^°\A,t) no longer has a power-law 
tail and instead takes on a simple gaussian form Eq. ( p.!9| ) sketched as the dashed line in 
Fig. 7. Comparing the solid and the dashed lines, we clearly see that the biggest effect of 
the disorder is to shift the distribution to the small-A end. From Fig. 7, it is clear that most 
of the samples have A < r^. In fact, the likely end point separations for typical samples 
is given by the peak of the distribution which occurs for A < a^. This suggests that most 
samples have a unique ground state (and a unique equilibrium state at low T), like the one 
sketched in Fig. 2(a). Only a small fraction of the samples, of the order 

U t ( T ) = r d d A W t (A,r) ~ rA(r/i), with 9 = 2( - 1, (3.23) 

J tC 

have almost degenerate ground states that are separated by a distance of O(r^) at a distance 
r 3> a T from the fixed end. These are the large scale active droplets excitations depicted 



in Fig. 2(b). It has been argued in Ref. [|L8] that ( > 1/2 for the low temperature phase 
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of the directed polymer in any finite dimensions. If this is indeed the case, then 9 > 
and from Eq. ( |3.23|) the large droplets are rarely encountered in the thermodynamic limit 
r — > oo. However, the rare occurrence of these droplets can dominate disorder- averaged 
thermodynamic quantities due to the diverging sizes of active excitations when they do occur. 
For instance, it is straightforward to verify that the mean square end point fluctuations, 
Ct(t, t) = |A 2 ~ r, as required by the statistical tilt symmetry Eq. ( |2.24j ), is dominatedby 
the large droplets in the tail of the distribution Wt(A, r), i.e., those with A ~ r^. 

At this point, we should reexamine the identification of Wt(A,r) with the droplet dis- 
tribution: In doing so, we made an implicit assumption that configurations with multiple 
almost-degenerate ground states give negligible contributions in the thermodynamic limit. 
This is clearly wrong in the absence of disorder. To check the assumption in the presence of 
disorder, we can compute the probability of having a triply-degenerate ground state, such 
as the one shown in Fig. 8. Extending the method described in the above sections, we can 
consider the probability that Pi(y, t \x, t) has three large peaks at the end point t = 0. 
(The same can be applied to degeneracies in the bulk.) This probability can be obtained 
from the correlation function P\P\Pi. The connected part of this correlation function of the 
distribution is formally given by G31 = 5 3 F/5J 3 and can be computed from the knowledge 
of G and Ti^ (see the Appendices). The resulting probability of yi, $2, and $3 all separated 
by 0(A) (Fig. 8), is of order A 2d W 2 (A,t) as one would naively expect from the droplet 
interpretation. Therefore, the occurrence of triply-degenerate ground states with three end 
points all 0(t^) away from each other is very rare, of 0(t~ 26 ). One can therefore neglect the 
effect of such multiply-degenerate ground states for most observables as long as 9 > 0, i.e., 
if C > 1/2. 

The analysis we have carried out yields information, as we have seen, about the proba- 
bility of finding large thermally active excitations, i.e., those with excitation free energy SF 
of order T or smaller. In contrast, the natural quantity that arises in the phenomenological 
scaling theory |L8] is the probability Pd(A, SF) of finding an excitation with displacement of 
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order A and free energy SF ~ A" where a is the exponent which controls the magnitude of 
the free energy variations as the end point (y, 0) is moved through distances of order A <C & 
(Sec. II. D). Since most of the path will change when A ~ ft, then t^ a ~ t 9 is the scale of 
variations of the free energy (when A > and also the variations from sample to sample. 
The basic scaling hypothesis is that all free energy scales on length scale r scale as r and 
all displacement A scale as r^, so that for <C A <C 

P D (A, 6F)d(6F) ~ ^p D (5F/A«). (3.24) 

From balancing the change in the mean free energy as a function of y (which is just 
K\y — x\ 2 /(2t) from the statistical tilt symmetry), with the random part of the free energy 
as y is varied, the exponent identity 

9 = 2C - 1 (3.25) 



is obtained. It was argued in Ref. that the natural form of the scaling function p£> is 
that it goes to a constant for small argument so that for T <C A", the probability of an 
excitation being thermally active is of order T/A a . Although "natural" and consistent with 
numerical studies, this does appear as an additional assumption. We now see, however, that 
the present calculation tells us that it is correct since we have found that 

P D (A,5F = T)~ W{A\t)d d A' ~ A~ 2+c , (3.26) 

J A 

which agrees with the conjectured form. The exponent 6 for active excitations with A ~ 
is therefore equal to 9. We have thus shown that the exponents describing the probability 
of the typical and the rare active excitations are the same! This result is rather important, 
since its analog has been assumed in a variety of random systems. 

Before discussing the physical response functions, we briefly consider what would happen 
if the wandering exponent ( were 1/2. It has been suggested, on the basis of results on the 
Calyley tree |Tj| , that this might occur for the pinned phase in sufficiently high dimensions. 



It also occurs at the critical point, separating the high and low temperature phases for 
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d > 2 ||32|| . If ( = 1/2, but the displacements at a distance r from the fixed end grow loga- 
rithmically faster than r 1//2 , then 9 would be zero but there will be logarithmic growth of the 
free energy scale with length scale and, presumably, a (logg)^ cusp in W(q, r) which would 
yield a distribution of the form W(A,t) ~ A~ d (log A) y in the scaling region. In this case, 
configurations involving multiply degenerate ground states could contribute substantially to 
disorder averages, and the droplet picture would no longer be simple. If ( is strictly 1/2 
with no logarithms, the behavior would be even more complicated. Exact calculations on 
the Cayley tree ]13| where this occurs shows that, indeed a finite fraction of samples have 
nearly degenerate ground states with very little overlap between them. However, many of 



the properties on the Cayley tree are pathological [|18[] because two paths never meet again 
once they branch. At the critical point of the polymer's glass transition for d > 2, the 



free energy scale is of order T because the transition occurs at a finite temperature [32 



Depending on the definition of 8, it may or may not include logarithms |Tj|. In any case, 
our considerations based on low free energy excitations must be modified as in Refs. [18| 
and PB|. The physics of this glass transition is very rich and will be left for future study. 



2. Anomalous Linear Response 

We now analyze the distribution of the linear response \. As explained in Sec. III. A, 
we obtain a power law tail for D t (x,r) from assuming that x is dominated by individual 
large active droplets of appropriate sizes. What were ignored were contributions to x horn 
multiply-degenerate ground state configurations such as the one sketched in Fig. 8. But we 
have shown in Sec.III.B.l that the probability of the triply degenerate configuration is of 
order t~ e lower than that of the doubly degenerate configurations. (Higher order degeneracy 
should be even rarer.) On the other hand, contributions to x from a given doubly or triply 
degenerate configuration are of the same order (~ A 2 /T). Thus it is reasonable to neglect 
the contribution to D t (x, t) from higher order degeneracies in the large r limit. The resulting 
distribution, which depends solely on the droplet excitations, should then be similar in form 
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to Wt(A ~ t) as sketched in Fig. 7 (solid line), with a regular part for small x an d a long 
tail described by Eq. ( |3.19|) for \ ^ ( a r/ ^) 2 - [Again, this would not hold at the critical point 
or in d = oo where ( = 1/2 and 6 = 0. In these situations, the multiply-degenerate ground 
states could contribute to % a ^ the same order, and the distribution function D t (x, T ) might 
be more complicated. It will be interesting to explore these special situations, for example 
by analyzing the directed polymer on a Cayley tree, or by studying the behavior at the glass 
transition.] 

The distribution D t iXi T )i an d in particular, the distribution of the linear response to a 
uniform field 

D( X ,t) = D t ( X ,t) (3.27) 

can be obtained experimentally, by following the domain wall of a planar random Ising 
magnet in the ordered phase, or by following the end points of flux lines in type-II supercon- 
ductors close to H c i. This may not seem practical at first glance, as many samples would 
seem to be required. However, thanks to the special statistical tilt symmetry present in this 



problem (see also Ref. []33|), the distribution D(x,t) can actually be obtained from a single 
sample. To see this, we just need to realize that the directed polymer only occupies a very 
limited volume of a given sample. To be more specific, the polymer of length t is essentially 
contained within a d + 1 dimensional "cone" (actually, a distorted paraboloid) of length t 
and radius p about the pinned end. So if we apply a uniform tilt field h [Eq. fl2.21|) with 



to = 0], the polymer will be biased, on average, to tilt by an angle = h/n. It is therefore 
forced out of the cone for 

h > 0(t f_1 ). (3.28) 

Once outside of the cone, the polymer encounters a random potential completely different 
from the one within the cone, since the random potential inside and outside of the cone are 
uncorrelated. Due to the statistical tilt symmetry, the random part of the new potential the 
polymer encounters is statistically the same as the one inside the cone. In other words, the 
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new potential is just like another realization of randomness or another sample except for the 
trivial mean tilt energy. Therefore, disorder averages can be obtained simply by monitoring 
the linear response \ as a function of the tilt angle which is varied by changing the applied 
field h over a range 5h ^> £* _1 . 

Recently, Mezard |2IJ reported such a numerical study of the polymer's response as 



a function of the applied tilt field h, for a fixed realization of the random potential in 
1+1 dimensions. Qualitatively, the position of the free end y(h) = (£(0)) h is found to be 
independent of h for a range of the applied field. But y(h) can "jump", by a distance of 
0(t 2//3 ) for a polymer of length £, if the applied field is increased by 0(t -1 / 3 ). As a result, 
the response x(h) = dy(h)/dh exhibits a series of sporadic sharp peaks, with strengths of 
0(£ 4//3 ), separated by distances of 0(£ -1 / 3 ), yielding x ~ t on average. (Here the average is 
performed over the applied field, i.e., by j- dh'x(h') for h ^> t -1 / 3 .) 

The above qualitative findings are of course consistent with the droplet picture described 
in this paper. As expected, a directed polymer in a typical realization of the random potential 
(or at some given tilt field) has an unique ground state. Typical excited states with distances 
0(t^) away are of O(t ) higher in free energy. This "energy gap" prohibits any large scale 
response to a very small increase in the field h. Since the applied field supplies an energy 
h ■ (£(£) — £(0)) ~ tit 1 *, the "energy gap" can be overcome for ht^ > t e . This is just the 
condition Eq. (|3.28 ) using the exponent relation Eq. (|3.25|) . In 1+1 dimensions, we have 



£ = 2/3. Thus the energy gap is overcome and the end point jumps when h ~ t -1 / 3 . By 
the statistical tilt symmetry, the process repeats itself (statistically) once the polymer's end 
point "jumps" by 0(t^). This then produces the sporadic jumps in y{h) [or peaks in xQ 1 )] 
observed in Ref. [pOfl . 

One might wonder why there are not many jumps of size A < t^, involving just a 
section of length A 1 ^ near the free end. These typically cost energy of order A a and gain 
field energy of order ShA. Thus naively, they will occur when 5h ~ l/A 1- ". However, 
since this increases with decreasing A, small jumps will almost alway be preempted by large 
jumps. They will only occur if the excitation energy is anomalously small [of order A/t 1 ^^ 
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which occurs with probability (A/t^) 1- "]. Therefore most of the end point motion as a 
function of h will occur in jumps of Oft). 

It should be possible to get some analytical information on the distribution of jumps from 
the methods of the paper, but the details are complicated. One would need information on 
joint distribution of end points with slightly different tilt fields h\ and h 2 , but the same 
random potential, such as 



0|f, t; h 1 )P 1 (y 2 , 0\x, t; h 2 ). (3.29) 

Conceptually, studying the response \ of a system by slowly varying an uniform bias h 
is similar to recent studies of self-organized criticality in cellular automaton models of sand- 
piles [0,0 and earthquakes []3~6|1 : Tilting the polymer corresponds to tilting the sandpile or 



increasing the loading force on a fault in earthquake models. Anomalously large responses 
(the sporadic peaks in xft)) exhibited by the polymers are like the "avalanches" encountered 
in the cellular automaton models. Of course, there are fundamental differences: The present 
study addresses equilibrium fluctuations of a disordered system at low temperatures, while 
the cellular automaton models mimic the nonequilibrium dynamics at zero temperature. In 
particular, if a directed polymer in a random potential with relaxational dynamics at zero 
temperature is driven by increasing the tilt bias adiabatically, the spectrum of avalanches, 
i.e., jumps, will be very different from the equilibrium behavior discussed above: There 
will be many small avalanches and only occasionally large ones, as analyzed for higher di- 
mensional driven manifolds by Narayan and Fisher ||37|| . In addition, as h is increased, the 
polymer will stay in configurations that for a while become more and more strained with 
the free end advancing much faster than the midpoint, in contrast to the equilibrium case. 
Nevertheless, the two problems share the key feature of a broad distribution of large rare 
"events". Intuition gained from studying the equilibrium behavior may therefore provide 
useful hints for understanding the nonequilibrium problem. 
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IV. DISCUSSION 



In this paper, we have presented a systematic investigation of some properties of the 
large scale, low energy excitations of directed polymers in random media. The large scale 
thermally active excitations occur in a small fraction of samples (or regions of a sample) 
which have almost degenerate ground states that differ by large displacement in the region. 
We showed the probability of occurrence of such samples, with excitations on scale t, is 
of order t~ e in the limit the polymer length t — > oo. Nevertheless, many of the disorder- 
averaged thermodynamic properties are dominated entirely by these rare samples. This is 
because in most of the samples, a polymer is essentially frozen in a unique ground state, with 
the largest scale excitations having energy of order t d hence not being thermally active. The 
rare samples with large active excitations make up for their small numbers by displaying 
anomalously large fluctuations which dominate the mean. We computed the distribution of 
the almost degenerate ground states for a polymer with one end fixed. While the distribution 
is peaked at small separations, we found a long power-law tail describing almost degenerate 
configurations that are far apart, with separation A at a distance r from the fixed end 
ranging from zero to the maximum important extent (with rapid decay for A ^> r^). This 
distribution of almost degenerate configurations also provides information on the distribution 
of excitations on scales much less than the system size which dominate many of the properties 
of large systems. 

The analysis was carried out by exploiting the known mapping of the directed poly- 
mer to the noisy Burgers equation, which converted singular, rare fluctuations to smooth, 
hydrodynamic fluctuations. We found that the interesting polymer distribution functions 
corresponded to various response functions of the Burgers' equation, and the tails of the 
distributions, which are generally difficult to obtain, are in this case just the hydrodynamic 
limits of the response functions and therefore much more straightforward to analyze. In the 
context of the Burgers' equation, the power law tails are merely the "long time tails" that are 
well known in hydrodynamic systems. Most of the difficulties with analyzing the properties 
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of the pinned phase of the directed polymer by conventional methods are due to its control by 
a non-trivial randomness-dominated zero-temperature fixed point at which temperature is 
"dangerously irrelevant" , meaning that interesting correlation functions cannot be obtained 
by setting T to zero. The main advantage of the approach via the noisy-Burgers' equation, 
is that the problem is converted into a much more conventional one in which fluctuations 
play the dominant role and field theoretic methods can be used. The physics of minimiza- 
tion of the coarse grained free energy to find "ground states" and then treatment of the 
anomalous rare fluctuations around the ground states miraculously appear as a consequence 
of the existence of a conventional fixed point of the hydrodynamic problem. 

Our findings support the main conjectures of the droplet theory |18| , |19[ . What we have 
demonstrated in this paper is a formal way of obtaining the statistics of the active droplets 
and various direct physical consequences of their distributions. As should be clear from the 
derivation of the distribution functions, a crucial ingredient which leads to power-law tails of 
distributions is the existence of the statistical tilt symmetry. It is thus tempting to attribute 
the scale-invariant distributions to some analog of the Goldstone modes associated with the 
tilt symmetry for the ensemble of disordered systems which is "spontaneously-broken" by the 
selection of the configuration for each specific sample. It would be important to understand 
whether in general such a symmetry is sufficient to ensure the existence of a power-law 
tail of large rare fluctuations in disordered systems controlled by randomness-dominated 
zero-temperature fixed points. 

The large rare thermally active excitations in the pinned phase of random directed poly- 
mers are closely analogous to the droplet excitations proposed to describe the low temper- 
ature phase of Ising spin glasses. Their existence in spin glasses was argued to depend, as 
in the present case, on the existence of a non-trivial zero temperature fixed point which 



is stable to thermal fluctuations ||10|| . This stability to fluctuations — and thus the whole 
picture — occur only when the free energy exponent 9 is positive (or perhaps zero with 
logarithms). But what is the analog of the statistical tilt symmetry for spin glasses ? In 
the Edwards-Anderson model, there is a statistical local gauge symmetry associated with 
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changing the sign of one spin and all the exchange couplings involving that spin. This was 
used, implicitly, to argue for the existence of large rare thermally active excitations in spin 
glasses JTIJ . It plays a role analogous to the statistical tilt symmetry of the directed polymer. 
Note however that the statistical gauge symmetry is not exact in more realistic models of 
spin glasses in which ferro and antiferromagnetic exchanges are not equally likely. Never- 
theless, it is believed that in the spin glass phase, the deviations from the exact symmetry 
are irrelevant on long length scales. A similar phenomenon occurs for directed polymers 
on a lattice in which there are short-distance correlations in the random potential along 
the polymer. In this case, the statistical tilt symmetry is not exact but numerical evidence 
support the conjecture that the deviations from it are irrelevant at long length scales. 

It is finally useful to compare this study with a recent uncontrolled field-theoretic anal- 
ysis using replicas. In principle, distribution functions such as W t (A,r) can be obtained 
using the replica method [|L7|| . However, difficulties arise in trying to take the n — > limit 
which is not uniquely defined from the positive integer n. Mezard and Parisi have pro- 
posed a gaussian variational ansatz |I7[ to treat the directed polymer problem. Within 



this framework, it is easy to see that replica-symmetry must be broken in order to recover 
the rare fluctuations. In such a treatment, the branching tree structure of ground states 
as a function of their end points in physical space (Fig. 6) is replaced by an abstract hi- 
erarchical "family" tree structure in replica space. The key to the power-law tails of the 
distribution functions is the overlap of "distant family members" , the u — > limit of the 



function [er](ii) in the notation of Ref. JT^]. Within the gaussian ansatz, a nontrivial tail 
is obtained if [a\{u — > 0) — > u 2 ^ e . Thus the replica-symmetry breaking scheme seems to 
be an approximate way of obtaining rare fluctuations within the replica formalism |39|| . As 
such, it may be a useful tool for other problems. However, the variational approximation is 
uncontrolled and the physics enters in a very roundabout manner. It is thus not at all clear 
why this approach has any advantages over a direct phenomenological scaling approach as in 
Ref. |T8[], renormalization group calculations which are exact on the hierarchical lattice [fTJ 



or approximate "functional" renormalization group studies |7,38|, all of which have the real 
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virtue that the physics appears directly. Furthermore, it is not clear if any real physics 
is associated with the notion of "replica-symmetry breaking" |39[]. As our study implies, 
the power-law tail of the rare fluctuations are a rather general consequence of the statistical 
symmetry of the problem; the ill-defined concept of "replica-symmetry breaking" is of course 



not needed at all in our analysis [3U| ! However, the simple structure of the free energy func- 
tional via the noisy-Burgers' equation, relies strongly on the one-dimensional nature of the 
directed polymer. It will be very interesting to see if our method can be extended, using the 
formalism of Appendix A, to study higher (but finite) dimensional systems such as oriented 
manifolds in random media 0. 

It may be comforting to some to find that the important physics of the ground state, 
- in particular power-law distributions of rare, anomalously large fluctuations — can be 
retrieved both by analysis in physical and replica space. For the directed polymer, at 
least, there are no substantial disagreement between these two approaches. It is possible 
that controversies regarding the nature of the low-temperature phase of finite-dimensional 
spin-glasses may have a similar resolution, with the many thermodynamic "states" of the 
infinite-range Sherington-Kirkpatrick model becoming, in finite dimensions, just the large 
scale, low free energy droplet excitations of the phenomenological scaling approach. 
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APPENDIX A: A FIELD THEORY FOR RANDOM SYSTEMS 



In this Appendix, we consider a general field theoretic description of manifolds in random 
media. In particular, we study a disordered system described by a Hamiltonian Ti.[£(z),r)] 
which is defined on the D-dimensional manifold z G K D , with £ e a (^-component "order 
parameter", and r](£,z) a quenched random variable (in our case the potential) distributed 
according to some distribution p[rj\. The one-point restricted partition function is 

Z(x , t_;r ] ) = Jv [e] 5 d (!*(£) - f) e^WT (A1) 

Various distribution functions can be obtained by adding a source term 

H^H + Jd D z J(i,z), (A2) 

and then differentiating the free energy, 

F(x,t_;J;r]) = —TlogZ(x, t; J; 77). (A3) 

For instance, 

SF(x, f ; J) 



^(faO-fi)), (A4) 



<5J(£i,£i) 

which is the disorder-averaged conditional probability of £(£1) = x*\ given = x, and 



--(#* (fao - ^) ^ (fob) - f 2 )) lc; , (as) 



5j(fi,ii)5j(f 2 ,i 2 ) T 

where the notations for thermal and disorder averages are the same as those used in Sec.ILA, 
and the superscript (c) denotes the connected (or truncated) correlations. We will generally 
be interested in the n th order distribution and correlation functions, 

Sj3(x i; . . . ; x„|x) = - 7 ^f (X; /y, V ( A6 ) 
oJ(xi) . . .dJ(x n ) 

and 

0g(x i; . . . ; x n ) = F( Xl )...F(x n ), (A7) 
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where we have used the shorthand Xj = (a?j, t-). Clearly, translational symmetry in space 
(afs) is restored upon disorder average. In addition, if we consider a very large manifold 
where the boundary of z is far from any of the distances — t , then we also have trans- 
lational symmetry in the internal coordinates. For a semi-infinite polymer, this leads to 
Oii(y,t — ti\x,t) = G(x — y, t — ti), which is just the one-point distribution function 
Pi(y, ti\x, t). Also, GqI(x, t; y, t) is related to the free-energy correlation function Cp{x—y, t) 
in Eq. (|2.28| ), and £?£i($i> *ij V2, ti\x, t) = G2,x(x — yi,x — y*2, t — t\) gives the droplet distri- 
bution through Eq. (|3.5f) . 

To carry out the average over disorder, it is advantageous to eliminate the random 
variables t](x) in favor of a probability density for the field F(x) itself, i.e., 

E[F, J] = J V[ V }5(F + TlogZ[J;r]})p[r]}, (A8) 

so that the mean free energy is 

F(x; J) — J T>[F)F(x)Z[F, J). (A9) 

It will be more convenient to write S in terms of its (functional) Fourier transform S[F, F] 



as 



E[F, J} = J V[iF] E[F, F]e/ rfxJ(x)FW . (A10) 

Note that Eq. ( |A10|) is similar to the dynamic generating functional of Martin, Siggia and 
Rose [|31],(|l[ with F being the "response field" . For a generic random system, the free energy 
functional E[F,F] will in general be rather complicated. However, the one-dimensional 
nature of the directed polymer can be exploited |||3| to derive a functional that is local in 
F and F. The explicit form of F] is shown in Appendix B. 

After adding a source term to Eq. (|A10 ) to generate correlations of F, we obtain the 



generating functional 



Z[J, J}= J V[iF]V[F] E[F, F] exp J dx (j(x)F(x) + J(x)F(x 



(All) 
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The distribution and correlation functions of interest can now be obtained from functional 
derivatives of $[J, J] = log Z[J, J], e.g., 



^(yi; • • -;ym|xi; . . . ;x„) 



SJ(yi) . . .5J(y m )5J(xi) . ..8J(x n ) 
F(y 1 )...F(y m )F(x 1 )...F(x ri ) (C) , 



where the disorder average (overbar) is now performed using Z[J = 0, J = 0]. 

Note that the normalization of J] requires Z[J = 0, J] = 1. Hence all of the 



moments G m o = 0. To obtain the other moments, we introduce the vertex functions [42 



8 m+n T[F,F] 

7m,n(yi; • • • ; y m |xi; . . . ; x n ) = — = = — , (A12) 

5F(y 1 )...8F(y m )8F(x 1 )...8F(x n ) 



where V[F, ~F\ is the Legendre transform of $: 



T[F,F] = J dx 



F(x)J(x) + F(x)J(x) 



$[J,J], (A13) 



with J and J eliminated from the above expression using 



— 5$ _ 5$ 

F = TT F = W (A14) 

It is a simple exercise to show that all of the correlation functions G$ n can be obtained from 



the vertex functions 7 mj „ The advantage of introducing the vertex functions is to allow 



the derivation of Ward identities from symmetries ||42|| . This will be done in Appendix C 



(c) (c) 

With g^' = 0, the expressions for Gmi are particularly simple: Only the connected parts 
survive (we will therefore drop the superscript (c)), and all of the higher order functions G m ,i 
can be obtained from the combination of 7i, m 's and Gi,i- For instance, 

J rfx 0i,i(y|x) 7 i,i(x|y') = 8(y - y'), (A15) 

thus 7i f i is the inverse "propagator" , and 

&,i(yi;y2|x) = - / rfx , rfy' 1 c/y2^i i i(x / |x)7 1>2 (x , |y , 1 ;y2)^i,i(yi|yi)^i i i(y 2 |y2). (A16) 
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This leads to Eq. ( |3.6|) for the directed polymer with (/^(x'lx) = G(x — x'), {?2,i(yi; Y2|x) = 
G 2 ,i(x -yx,x- y 2 , r), and 7i, 2 (x'|yi; y 2 ) = r li2 (f - y[, f - ^; x' - j? 2 , t' - t' 2 ) due to the 
translational symmetry in a semi-infinite polymer after averaging over the disorder. The 
result Eq. ( |A16| ) is conveniently represented by a tree diagram (Fig. 6). A higher order 
distribution function Q m ^ typically contains a tree diagram involving only 7 12 and 
plus diagrams involving higher order vertex functions 7i, m ', with 2 < m! < m. An even 
larger class of vertex functions are needed to generate the correlation functions, {?o,n's. The 
description of these is more complicated, and we will not discuss them except to note that 
for the directed polymer in 1 + 1 dimensions, there exists identities (fluctuation-dissipation 



theorems) linking Q n and Q n ^i [p6| , [43" 
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APPENDIX B: MAPPING TO STOCHASTIC DYNAMICS 



In this appendix, we describe a mapping of the random directed polymer to a stochastic 
dynamics problem — the noisy Burgers' equation — first noted in Ref. @- We show how 
this mapping can be used to extract information about the directed polymer, and thereby 
provide a concrete example of the general field theoretic method described in Appendix A. 

Let us consider the doubly-restricted partition function 

Z 2 (x , t; x>, 0) = J D[£] 6 d (i(t) - x) 6 d (f(0) - x') e~ n / T (Bl) 

where (x, t) and (x 1 , 0) are the end point locations of a directed polymer of length t. We 
see that Z 2 is the Boltzmann weight of "propagating" the polymer from (x',0) to (x,t); it 
satisfies the "diffusion" equation 

T^-Z 2 (x, t; x', 0) = %-%Z 2 - v{x, t)Z 2 , (B2) 
at Ik 

where t is interpreted as "time" and (af, 0) is merely a parameter, which enters through 
the initial condition Z 2 (x,t = 0;x',0) = 5 d (x — x'). Since Eq. (|B2| ) is linear, we can 
easily integrate over the coordinate x' to obtain the same diffusion equation for Z(x,t) 
I d d xZ 2 (x, t; x', 0), but now with the initial condition Z(x, 0) = 1. Then the free energy of the 
one-point restricted polymer, F(x,t) = — T log Z(x,t) satisfies the well-known noise-driven 
Burgers' equation HiTj[l2| , 

ap T - 1 

with the initial condition F(x, 0) = 0. Note that the path integral description of the directed 
polymer, i.e., Eqs. (|2.1| ) and (|2.3| ), is very different in form, from the dynamic description 
Eq. ( p3|) |12| . The task of global optimization presented in the original directed polymer 
problem is turned into a step-by-step recursive process: Eq. (|B3| ) relates the free energy of 
a polymer of length t to one of length t + dt and is thus a transfer matrix description of the 
problem [[4(] . 
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The noisy-Burgers equation Eq. ( p3|) itself has been the focus of several recent stud- 
ies ||43|| . However, despite much effort, no systematic solution has been found. The lone 
exception is in 1+1 dimensions, where Eq. (|B"3"|) can be thought of as the d — 1 limit of an 



anisotropic driven diffusion equation p4|.|45| 



2k Ik 

with (rj'(r, t)rj'(^, t')) = -2DV 2 5 d {f- 7)5(t - t'), (B4) 



where x is one particular component of f. Eq. ([BJ) is equivalent to Eq. (|B3|) in d = 1 
with the identification = dF/dx and rf = dr]/dx. [See Ref. [[16|] for a detailed discussion 
of the connections between Eq. (|B3|) and Eq. ([B4|).1 Eq. (|B4j) has also been investigated 



in detail [|1[ , |5| . It can be solved perturbatively via a d = 2 — e renormalization group 
expansion. Furthermore, all of the scaling exponents can be obtained exactly in all d, owing 
to the existence of a "fluctuation-dissipation theorem" which we shall not elaborate here. 
One obtains the directed polymer results £ = 2/3 and 6 = 1/3 which we have referred to 
several times in the text. Not much is known analytically of the noisy-Burgers' equation 
away from 1+1 dimensions. Recently, an uncontrolled mode-coupling approach [f47|f48l along 



the line of Ref. []26[ has been applied to solve Eq. ( |B3|) in 2+1 dimensions, yielding exponent 
values that are rather close to those obtained from numerical simulations |43[]. However, a 
controlled solution is yet to be developed. 

Although quantitative solution (e.g., the exponent values) of the noisy-Burgers' equation 
are not available for d ^ 1, we can still learn much about the qualitative behavior of the 
directed polymer, e.g., the rare fluctuations, by exploiting this mapping to hydrodynamics, 
and especially by exploiting symmetries, as we now describe. From the standard dynamic 



field theory formulation of Langevin dynamics upon averaging over r] f3^||l| , one can turn 
Eq. ( P33|) into the free-energy functional used in Eq. (|A11|) , 



Z[F,F} = e- s ^ F \ (B5) 

where 
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S[F,F] 



dtd d x 



DF 2 (x,t) + F(x,t) ld t F(x,t) - ^V 2 F+ ^-(VF) 



2k 



2k 



(B6) 



Since the disorder potential rj has already been averaged out, translational symmetries in 
both x and t are restored in the limit of large t, corresponding to the hydrodynamic steady 
state. It is then convenient to go to Fourier space, where 



S[F,F}- 



du d d k 



+ 



/ 



2tt (27r) d 

duji du>2 d d k\ d d k 2 



D\F(k, u)\ 2 + ( —k 2 - iuj ) F(-k, -u)F(k, uj) 
2k. 



(B7) 



2tt 2tt (27r) d (2ir) d V 2k 

The functional S[F,F] is called the dynamic functional in the study of stochastic dynam- 
|T[. In the absence of randomness, i.e., for D = 0, this is just the "bare" generating 



ICS 



functional of the vertex functions, 



S[F, F] = r®[F,F]. 



(B8) 



From the definition of the vertex functions (|A12p , we have 

1 



r[F,F] = £ 



m\n\ J(k,u>) 



T-i pm pn 



(B9) 



where T mtn is the Fourier transform in x, y and t of the function 7 m>n in Eq. (A12). Com 
paring the above to Eq. (|B7]), we can easily read off the bare vertex functions, 

?(°) . T 



r i,i(A;,a;) = —k -iu, 
2k 

»(o) r _ 1^ ^ 
r i2 (/ci,a;i; £2,^2) = &i • h, 

K 

^(0) 

r 2,o( fc ^) = °- 



(BIO) 
(Bll) 
(B12) 



The relevance of the random potential in the pinned phase makes the renormalized func- 
tional r[F,F] much more complicated. However, we will show in Appendix C that, as a 



consequence of the statistical tilt symmetry, the scaling behavior obtained from rf 2 and 
are the same. This permits us to use the bare vertex Tf 2 in place of the full vertex Ti 2 to 
compute the tail of the droplet distribution W(A, r) in Sec. III. A. 
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As shown in Sec. III. A, W(A,r) is obtained by adding a source term J to the random 
potential i], and then following the changes in F. In the context of stochastic dynamics, 
this procedure amounts to probing the "response" of the hydrodynamic field F to a small 
"perturbation" J added to the right-hand side of Eq. ( |B3l) . We see that the distribution 
functions 5F/5J = G and 5 2 F/5J 2 = G 2 . i are simply the linear and nonlinear response 
functions of the noisy- Burgers' equation. Furthermore, the tails of the distribution functions 
are similarly related to the hydrodynamic limits of the response functions. Thus the mapping 
to hydrodynamics provides us with a convenient way of approaching the tails of distributions. 
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APPENDIX C: STATISTICAL TILT SYMMETRY AND WARD IDENTITIES 



In this Appendix, we review the statistical tilt symmetry which was used in Ref. |2P 



to prove the result \ = t/k. The analogous symmetry in the noisy-Burgers' equation is 
the Galilean invariance. It has been used in Refs. [0,0 to derive the exponent identity 
6 = 2£ — 1. Here we shall derive the consequences of this symmetry, i.e., Ward identities, 
more formally. The result is used to compute the droplet distribution in Sec. III. A. 

We start with the Hamiltonian 7~C[£,rj\ given in Eq. fl2.1|) . Consider the effect of an 
arbitrary bias of the form 

n h [l v] = H[l V] - [ dz h(z) ■ (CI) 

By a simple coordinate transformation £'(z) = £(z) — t(z)/k, where f(z) = J z dz' h(z'), we 
obtain 

H h [l rj\ = H[t V\ ~ ^ f dz' P(z'), (C2) 
2k Jo 

with a transformed random potential 

rj(£', z )=r][£'{z)+r(z)/K, z\. (C3) 

Similarly the free energy of the polymer with a point fixed at = x is given by 

F h (x, t; rj) = F(x - f(t)/K, t;rj) - ^- / dz h 2 {z). (C4) 

2k Jo 

Note that since the random potential 77 is gaussian and uncorrelated, the distribution of rj 
is also gaussian, characterized by its covariance 



fj(xi, ti)rj(x 2 , t 2 )= rj(xi, ti)T](x 2 , t 2 ) 

= 2D5 d {x 1 ~x 2 )5{t 1 -t 2 ). (C5) 

Thus disorder averages of a functional of F over rj and rj are identical. In particular, if we 
choose h(z) = h for t < z < t, and h(z) = for < z < t , then Eq. ( |CTD reduces to 
Eq. fl2.21|) , and Eqs. (|C4|) and (|C5|) imply that the disorder-averaged linear response is just 
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X--^F h (x,t) 



1 ~ to (C6) 



This result was first proved for a class of disordered systems in Ref. [f29| . 

In the following, we shall exploit this statistical symmetry more formally by considering 
in detail the case to = 0, corresponding to applying a uniform tilt field h. In this case, 
Eqs. ( |C3| ) and ( JJty become 



rj(£',z)=r][Z'{z) + hz/K, z), (C7) 

and 

h 2 

F h {x, t; rj) = F(x - ht/K, t; rj) - —t. (C8) 



From Eq. it is clear that no longer satisfies the noisy- Burgers' equation (|B3| ) 



However, the combination F^x, t) —Th ■ x does. This is because Eq. ( |B3| ) is invariant under 
the "Galilean" transformation 

F'(x, t) = F(x + vt, t) — kv ■ x, (C9) 

and 

rj'(x,t) =rj(x + vt,t), (CIO) 

as can be verified straightforwardly for v = h/n. We now use this Galilean invariance to 
derive some useful Ward identities. 

The statistical tilt symmetry implies that the generating functional of the vertex function 
T[F, F] is invariant under the Galilean transformation of the disorder-averaged fields F and 
F. In Fourier space, the transformation reads, 

F'(q, t{)= e^F(q, t) + ikv ■ V^ d (g), (Cll) 
F{qM)=^F{qM)- (C12) 

Thus we have the identity 
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6T 

Sv 



= / d d qdt x { i qhFi&tt) - K5 d (q)V^ 



5Y 



iqUF^tx) — 



ST 



(C13) 



where Eqs. ( pll ) and (|C12|) have been used to obtain the variations of F and F with respect 
to v. Taking functional derivatives of Eq. (|C13| ) with respect to F(k,t) and F(k — q,t 2 ), 
then taking the limit F, F — > 0, we obtain the following Ward identity in Fourier space, 



kV^J dt x Y^ 2 {q,t- t x \ k - q,t - * 2 )|f-»o = % - t 2 )f ltl (k,t - t 2 ), 



where the definition 



(C14) 



Ti, n (,qi, t -t 1 -----q n ,t- t n ){2<K) d 5 d {q\ + ■ ■ ■ + q n + k) 



= / d d xd d X\ ■ ■ ■ d d X Hq n £n 



S^+n) T [F,F] 



5F(x, t)5F(x 1, tt) ■ ■ ■ 5F(x n , t Tl 
is used. It is more convenient to write Eq. ( |U14 ) in frequency space. With 



ri tn (q\,ui; q n ,u n )27T5(uJi H hw„ + w) 

dtdh ■ ■ ■ dt n e~ ilJt - i ^-~ i ^ tn f l>n {q 1 ,t -tu---;q n ,t-t n ), 



Eq. ( |C14| ) becomes 







«V$Ti,2(g',/i; k-q,u- /x)|f >M -»o = -ik—T hl (k,u). 



From Eq. (|A15 ), we have 



--1 



r 1;1 (k,u) = G (k,u), 



(C15) 



(C16) 



where 



G{k,u) = / dtg(qf) 



In the simplest case where there is no disorder, we have 

s(o) T 

T 11 (k, uj) = -^-k 2 — iu. 



(C17) 



(C18) 



The Ward identity then reads 
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«V f r i 2 (g,//; k-q,u- A*)|«/*-rf) = (C19) 



yielding a "bare" three-point vertex 



x(o) _ 1^ ^ 

r 12 (fci,u;i; ^2,^2) = fei-fe, (C20) 



since Ti^ must be symmetric in k\ and fo. This vertex is exactly Eq. ( [BHD , which has been 
obtained directly from the dynamic functional S[F,F] in Appendix B. 

For the random problem at hand, the scaling form Eq. (|2.9|) or (|3.9| ) for G implies that 
the two-point vertex must have the following form: 

f 1>1 (k 1 u) = k 1 ^f 1>1 (u/k 1 ^) 1 ( C2i) 

where 

f 1A (s) = v 1 (s)-tsv 2 (s), (C22) 

with Ux, v-i being real, dimensionless scaling functions. Since g is real, v\ and v 2 must be 
even functions. Furthermore, / °° dtg{qt C: ) ~ q^ 1 ^ since g is rapidly decreasing for large 
argument. It then follows that 

lim^x(s) = const and lims ^2( s ) = 0. (C23) 

s^O s— >0 



The Ward identity ( |C15|) now leads to the relation 



^ — * 1 — # 1 //■ 

ri )2 (?,A*; k-q,u- yu)|g>^ = q ■ {k - q)v{u/k (C24) 

K 

where u(s) = i-j^Tx^s) is another dimensionless (but complex) scaling function, 

/\ /\ dvo .dv-\ /^,„_x 

KsH^«+«T- ! T' ( C25 ) 
ds ds 

The above limiting behavior of r 12 is as much as the Ward identity will directly yield. 
However, it does suggest a natural scaling form for the full spatial and temporal dependence 
of the vertex, 
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fi, 2 (^i,^i;^2,^2) = — £1 • k 2 V(u 2 /kl /c ,k 1 /k2,uJi/u;2), (C26) 

K 

where V is another complex, dimensionless scaling function, with V(s, 0,0) = v{s). 



The scaling form Eq. ( |C26| ) of the full vertex r^, together with mild convergence con- 
ditions on V, imply that the scaling properties obtained from the full vertex is the same 
as that obtained from the bare vertex (|C20| ). We now illustrate this by computing the tail 
of the droplet distribution, or equivalently the singularity in W(q,r — > oo), using the full 
vertex. 

In the limit r — > oo, it is convenient to write Eq. ([3 .10 ) in frequency space, 

W'{q, t -> oo)= — T / — r 1>2 (5, u; -q, -u)G(q, uj)G{-q, -u) 

= -- ? 2 -VC /; ( C27 ) 

K 

with 

-oo27TZ/ 1 2 ( S ) + S V 2 2 ( S )' 1 ' 

where we used Eqs. QC22Q and ( C26 ). The scaling form Eq. (|3.15| ) is valid as long as the 



integral I converges. In Sec. III. A, we used the "bare" T's corresponding to V — v\ — v 2 — 1. 
This gives the result 7 = 1/2. In general, we see that, given the properties of v\ and v 2 in 
Eq. ( |C23| ), the integral converges if 

V(s->0, 1,1) < s- 1 and V(s -> oo, 1, 1) < s u$(s). (C29) 

Unfortunately, V(s, 1,1) is not directly obtainable from the Ward identity, which only gives 
V(s,0, 0). However, we do not expect V(s,ki/k2,uji/uj2) to depend strongly on the ratio 
ki/k 2 and lo\/l02 in the region < ki/k 2 < 1 and < uji/uj 2 < 1. If we use V(s, 0, 0) = u(s) 
for V(s, 1, 1) with i/(s) given by Eq. ( |C25| ) , then the convergence conditions (|C29|) are 
always satisfied given the properties of V\ and v 2 in Eq. ( |C23|) . It is thus very reasonable to 



expect that (|C29| ) are also satisfied by V(s, 1, 1) itself, which then leads to the scaling result 
Eq. (ETI5D for W'{q, r), with w(s -»• oo) = I. 
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Finally we note that although the scaling form ( |C26|) for the vertex and the convergence 
conditions ( |C29| ) have not been derived from first principles, they are very plausible if an 
underlying renormalization group exists. This is especially true in 1+1 dimensions, where the 
strong coupling fixed point of the noisy-Burgers' equation can be described perturbatively via 
a d = 2 — e renormalizat ion-group expansion of the anisotropic driven diffusion equation ( |B~I| ) 
in d+ 1 dimensions as explained in Appendix B. The vertex function can be computed 
explicitly in the 2 — e expansion, and is indeed described by the scaling form Eq. (|C26| ). 
Away from 1+1 dimensions for the noisy-Burgers' equation, a perturbative renormalization 
scheme is not known. However none of the known results give any indication of the break 
down of simple scaling. Thus we expect that the hypothesized scaling form ( |C26| ) for r 12 is 
valid in all dimensions for the Burgers' equation and therefore also for the directed polymer. 
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FIGURES 
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FIG. 1. (a) Schematic of the ground state £*(z) of a typical directed polymer with one end 
fixed at (x, t) (the solid circle). The position of a segment at z = to (open circle) is at £*(*o) = y*. 
(b) The one-point distribution function Pi(y, to\x, t) as a function of y. It is sharply peaked near 
some sample-specific y* with \y* — x\ ~ r^, where r = t — to- The dashed line indicates the 
disorder- averaged distribution function Gt{x — y,r) = Pi(y,to\x,t), which is centered about x, 
with a width r^. This gives the distribution of y* from sample-to-sample. 
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FIG. 2. (a) Behavior with a configuration of the random potential that gives rise to two 
almost-degenerate "ground states" ££(2) and H(z) for a directed polymer with one end fixed at 
(x,t). The positions of the free ends of the two states are £a(0) = y* a and £t(0) = yl- (b) The 
end-point distribution function P\(y, 0|x, t) for the sample in (a) displays two peaks at y^ and y^. 
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FIG. 3. The dotted lines are the degenerate ground states Q( z ) an d Cb( z ) shown in Fig. 2(a). 
The solid line is the optimal path with both ends fixed: £(i) = x and £(0) = y, where y is between 
y* a and y£, which are separated by a distance of order A. 
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FIG. 4. The solid line illustrates a possible degenerate configuration of size r in the bulk for 
the optimal path shown in Fig. 3. Some energy barrier separates the two degenerate pieces of the 
optimal paths to y. These will affect the energy barrier to move the polymer from between its 
almost-degenerate ground states ending at y* and y\. 
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FIG. 5. Two polymers with ends fixed at a distance A apart typically merge at a distance A 1 ^ 
from the fixed ends. 
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FIG. 7. The distribution function of the active droplet of size A, W(A, r) (the solid line), is 
composed of two parts: A smooth piece for A < a^(T), the small scale cutoff, and a long power-law 
tail extending to "K. The dashed line is the gaussian distribution W^°\A,t) which describes the 
pure system (and the high temperature phase of the disordered system). It has a width of 0(r 1 / 2 ). 
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FIG. 8. A triply degenerate ground state configuration for a directed polymer in 2+1 dimen- 
sions, with the distance between each pair of end points being of order for a polymer of length 
t. 
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